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Abstract 

The one-loop contributions to the branching ratios for leptonic 
r decays are calculated in the CP conserving 2HDM(II). We found 
that these one-loop contributions, involving both neutral and charged 
Higgs bosons, dominate over the tree-level exchange, the latter one 
being totally negligible for the decay into e. The analysis is focused on 
large tan (5 enhanced contributions to the considered branching ratios. 
We derive a simple analytical expression for the one-loop contribution 
which holds in this case. We show that the leptonic branching ratios 
of r are complementary to the Higgsstrahlung processes for h{H) and 
have a large potential on constraining parameters of the model. In 
this work we provide upper limits on Yukawa couplings for both light 
h and light A scenarios and we derive new lower limit on mass of 
M H ± as a function of tan /3, which differs significantly from what was 
considered as standard constraint based on the tree-level exchange 
only. Interestingly we obtain also an upper limit on M H ± . For a SM- 
like h scenario, with heavy and degenerate additional Higgs bosons, 
one-loop corrections disappear. 
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1 Introduction 



The mechanism of electroweak symmetry breaking is the most important 
ingredient in the description of elementary particle physics. The Standard 
Model (SM) incorporates the Higgs mechanism that breaks the electroweak 
symmetry spontaneously through a neutral scalar field with non-zero vacuum 
expectation value. In the minimal version of this mechanism one scalar 
SU(2)l doublet is required, providing one physical particle: the Higgs boson. 
The search of this particle is one of the main aims of high energy physics 
and current searches at LEP exclude SM Higgs bosons with masses below 
114.1 GeV at 95% C.L. [1]. In this context, valuable information about the 
Higgs mass will come from analysis of precise measurements of electroweak 
observables. The result of these indirect searches gives an upper bound on 
the Higgs mass M HsM < 219 GeV at 95% C.L. [1], that is of great importance 
for future searches. 

Models with Two Higgs Doublets (2HDM) are the minimal extensions of 
the SM Higgs sector describing all high energy experimental data and pro- 
viding new and rich phenomenology. These models can also be interpreted 
as effective theories describing low-energy physics in models with beyond the 
SM physics at higher scale. This maybe the case of the Minimal Supersym- 
metric Standard Model (MSSM) with heavy supersymmetric particles. A 
CP-conserving 2HDM contains 5 physical Higgs bosons, two neutral scalars, 
h and H, one pseudoscalar A, and two charged Higgs bosons, H ± (see e.g. [2]). 
The LEP direct searches for these Higgs bosons are more complicated than 
in the SM due to the number of free parameters involved and, in particular, 
the existence in 2HDM of one very light Higgs boson can not be excluded [1]. 
On the other hand, LEP data excludes, for example, neutral Higgs bosons 
with masses below 40 GeV in the regime of large tan f3, equal 60 or larger [3]. 

In this context, indirect searches for 2HDM effects in electroweak observ- 
ables provide important information about the masses and mixing angles in 
the Higgs sector. For example, concerning the charged Higgs, a lower bound 
M H ± > 490 GeV can be set using indirect effects in b — ^ 57 [4], to be com- 
pared with M H ± > 75.5 GeV coming from the direct LEP searches [5]. In 
order to explore the whole parameter space, global fits using different elec- 
troweak observables p, Rb and b — > S7 [6,7] (and also (g — 2) M in [7]), have 
been made, constraining large regions of the parameter space and therefore 
giving valuable information for future searches. 

In this work, a complete study of one-loop 2HDM effects in the leptonic r 
decays is performed for large tan (3 and arbitrary Higgs spectrum, extending 
the results from [8] and [9]. 

It will be seen that this radiative effects in the branching ratios for r — > 
ez/ e z/ T and r — > fiu^Ur are larger that the 2HDM tree-level effects in the 
relevant regions of parameter space and experimental data will be used to 
derive new constraints for Higgs masses and mixing angles. 

The paper is organised as follows. The Sect. 2 contains a short description 
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of the 2HDM properties and of results of the experimental searches on Higgs 
bosons. In Sect. 3 the leptonic r decay data are compared with the SM 
prediction and the 95 %CL bounds for 2HDM contributions to the branching 
ratios are derived. In Sect. 4 the 2HDM contributions are parameterised. 
The one- loop 2HDM effects are computed in Sect. 5 while their numerical 
analysis is performed in Sect. 6. Finally, in Sect. 7 we derived the constrains 
on the 2HDM parameters coming from leptonic r decay data analysis and 
our conclusions are summarised in Sect. 8. 

2 CP conserving 2HDM Model II 

2.1 General properties 

The Two-Higgs-Doublet Model is the simplest extension of the Standard 
Model with one extra scalar doublet. It contains three neutral and two 
charged Higgs bosons. Here we consider a simple CP conserving version 
with a soft Z 2 -violation, assuming the Yukawa interactions according to the 
Model II, 2HDM(II), as in MSSM. In this model one of the Higgs scalar 
doublet couples to the up-components of isodoublets while the second one 
to the down-components. In this case there are 7 parameters describing the 
Higgs Lagrangian: four masses for h, H, A and H ± , two mixing angles a and 
(3 (used in form sin(/3 — a) and tan/3 = v 2 /vi), and the //-parameter, related 
to the soft-Z 2 violating mass term in the Lagrangian (u = 3?m^ 2 /2t>!t> 2 ). 
This //-parameter describes the Higgs selfcouplings if they are expressed in 
terms of masses. We stress that none of these selfcouplings are involved in 
this analysis directly. However our results are sensitive indirectly to the v 
parameter as this parameter governs the decoupling properties of the model. 

There is the attractive possibility of having a neutral Higgs boson h sim- 
ilar to the SM one, and all other Higgs bosons much heavier. This scenario 
can be realised in two ways, depending on the value of //-parameter. For 
large v the additional Higgs boson masses can be very large and almost de- 
generate, since all of such masses arise from one large parameter - v. It 
is well known that in such case there is decoupling of these heavy bosons 
from known particles, i.e. effects of these additional Higgs bosons disappear 
if their masses tend to infinity, eg. in the 77/1 coupling. At small v the 
large masses of such additional Higgs bosons arise from large quartic selfcou- 
plings (A) in the Lagrangian. Since these couplings are bounded from above 
by the unitarity constraints, so are the heavy Higgs-boson masses. Accord- 
ing to these bounds heavy Higgs bosons have to be, typically, lighter than 
600 GeV [10]. Therefore, in this scenario the additional Higgs bosons can be 
heavy enough to avoid direct observation even at next generation of colliders, 
although some relevant effects can appear in the interaction of the lightest 
Higgs boson (non-decoupling) [11-14]. 

Another interesting scenarios that will be intensively studied in this work 
are the ones with mass of h or A below the SM Higgs boson mass limit, 
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114 GeV. In particular, the light A scenario is especially relevant for the 
description of (g — 2) M data [15-17]. These scenarios are possible without 
conflict with the existing data within a 2HDM(II), since this model allows 
for low production rates for very light Higgs particles, as will be discussed 
below. 

The 2HDM (II) model is characterised by the couplings of Higgs bosons to 
the fermions and to the EW gauge bosons. Their ratios to the corresponding 
couplings in the SM, \) — 9j/9j M i are presented in Table 1. Note that for 
couplings to the EW gauge bosons V, 

(X h v) 2 + (X?) 2 + {xtf = 1, 

and similarly for the couplings to fermions [11]. Note also, that for each 
neutral Higgs boson (i) we have 

(xi + X d )Xv = 1 + X u X d - 
Table 1: Relative couplings, Xj = d'j/dj^ 1 in 2HDM (II) 





h 


H 


A 


Xv 

Xu 

Xd 


sin(/3 — a) 

Xv + ™t(3xv 
Xv -t&nPxv 


cos(/3 — a) 

Xv-^t(3 X h v 
Xy+tan/?x{/ 




— i cot (3 
—i tan p 


XW-H+ 


cos(/3 — a) 


sin(/3 — a) 






Note, that for large tan/? the couplings to the charged leptons (equal 
to the couplings to the down- type quarks Xd)> relevant for our analysis, are 
enhanced. 

In the last row of Table [T] the W ± H T 4>° couplings, with 0° = h, H, A, 
being of interest of this work, are presented. Here the ratios of such couplings 
to the SM Higgs boson coupling to gauge boson, x\v-h+ = 9w ± h^ 1 1 9w ' ■> are 
shown. 

It is important to notice the complementarity between the x\/ on one 
hand and x\v-h+ ( anc ^ Xd a ^ iar g e tan/5) on the other. 

2.2 Experimental constraints on 2HDM 

The most important constraints on the 2HDM(II) parameter space come from 
LEP direct searches for Higgs bosons. Concerning light neutral Higgs bosons 
production, there are three main processes within the energy range covered 
by LEP, namely, the Higgsstrahlung, e + e~ — » Z* — > Zh, the associated 
production, e + e~ — > Z* — >• hA, and the Yukawa processes, e + e~ — > ff — > 
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ffh(A). The two first processes are highly complementary, due to their 
dependence on (/3 — a), 

cr(e + e~ — > Z* — > Z7i) = sin 2 (/3 — a)<7sM(e + e~ — > Z* — > ZH$m) 

a(e + e- -> -> = cos 2 (/3 - a)a 5M (e + e~ -> -> ZH SM )\ (1) 

a{e + e-^Z*^ffh) = ( x h d ) 2 a SM (e + e-^Z*^ffH SM ) (2) 

a ( e + e -^ Z *^ffA) = ( x ^a SM (e + e-^Z*^ffHs M ), (3) 

where A = A^/[A^(12Ml/s + A zfc )], with A^- = (1 - m?/s + mf) 2 - 
4m 2 m 2 /s 2 , being the two-particle phase-space factor. 

The search for a light h through the Higgsstrahlung process, under an 
assumption that the light Higgs boson decays into hadronic states, has been 
performed in [18]. The results of this analysis set an upper limit on the 
product of the cross section and the corresponding branching ratio. It can 
be translated into an upper limit on sin 2 (/3 — a) as a function of Mh, shown in 
Fig. U (left) [18]. Therefore, the results of this analysis are compatible with 
a light h scenario (with mass below 114 GeV) if sin 2 (/3 — a) is small enough. 

Also upper limits on the cross section of the associated hA production 
process have been derived assuming 100% decays into hadrons [19]. These 
results can be translated into forbidden regions in the 2HDM(II) parameter 
space. In particular, these results highly constrain a scenario with both h and 
A light (the light Ahh scenario). In Fig. ^ the excluded (M^Ma) regions 
have been plotted [19]. A particular point is excluded in Fig. [T] (right) if 
it is excluded for 0.4 < tan/? < 40 (darker grey region), 0.4 < tan/? < 1 
(lighter grey region), and 1 < tan/3 < 40 (hatched region) for all values of 
a = ±7r/2, ±7r/4, 0. It is noticeable that a scenario with light h ( A) is not 
excluded if Ma ( Mh) is large enough. In particular, if sin 2 (/3 — a) = 0, 
LEP measurements are sensitive to this associated production if Mh + Ma < 
130 - 140 GeV. 

Finally, the search for a light Higgs boson has been performed through 
the analysis of Yukawa processes assuming that Higgs boson decays into r, 
if 2m T < Mh, Ma < 2rrib, or into 6-quarks, if Mh, Ma > 2rrib [3]. One of the 
results of this analysis is that M^a — 40 GeV are excluded for high tan/3 
(tan/3 > 60). We will discuss existing constraints together with new ones 
coming from our analysis in Sect. 7. 

Concerning the charged Higgs boson, direct searches at LEP through the 
process e + e" — > H + H~ have been performed assuming BR(H~ — > qq) + 
BR(H- -> tv t ) = 1. The lower bound M H ± > 75.5 GeV at 95% CL [5] 
was obtained. The Tevatron data set limits on mass of the charged Higgs 
boson as a function on tan /3, they are presented together with LEP results 
on Fig. |21 

Much stronger constraints on Mh± come from the indirect charged Higgs 
boson effects in b — > sj processes, if interpreted in 2HDM(II). This leads to 
a lower mass limit of 490 GeV at 95% for tan (3 > 2 [4] 4 . 

4 Recent analysis on B — > X s ^ predicts larger theoretical errors in the SM prediction 
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m^GeV/c 2 ) m h [GeV] 

Figure 1: Left: Upper limit on sin 2 (/3 — a) as a function of [18]; 
Right: Excluded (M^Ma) regions in the different ranges of tan /3 for a = 
±tt/2, ±tt/4, by OPAL [19]. 



In this context, important information on the available 2HDM(II) param- 
eter space is coming indirectly from the low energy precise measurements. In 
particular, from the Upsilon decay into h(A)j and g — 2 data, see eg. [17]. 
Also global fits have been performed, combining the results coming from the 
different electroweak observables p, Rb and b — > 57 [6, 7] (and also (g — 2) M 
in [7]), constraining large regions of the parameter space. Here, indirect con- 
straints of 2HDM(II) will be obtained using leptonic r decays data. The 
obtained results will be compared with direct search analysis coming from 
LEP and some low energy experiments. The implementation of leptonic r 
decay data in global fits will be performed elsewhere. 



3 Leptonic r decays: data versus SM predic- 
tions 

We consider the partial decay widths and branching ratios for the two leptonic 
decay channels of the r-lepton, namely 

t — ► ev e v T andr — > pPfjy T . (4) 

We will denote the corresponding quantities using superscript I, I — e and p, 
for example for the branching ratio we use Br 1 = Br(r — > WiV T ). 

The '04 world averaged data for the leptonic r decay modes and r lifetime 

and therefore a more conservative lower bound M H ± > 200 GeV [20] . 
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are [1] 

Br e \ exp = (17.84 ± 0.06)%, Br% xp = (17.37 ± 0.06)% (5) 
t t = (290.6 ±1.1) x 10~ 15 s. (6) 

Note that the relative errors of the above measured quantities are of the 
0.34-0.38 %, the biggest being for the lifetime. 

The SM prediction for these branching ratios can be defined as the ratios 
of the SM predicted decay widths to the total width as measured in the 
lifetime experiments, namely Br l \ S M = ^Ism/^xp — ^ 1 \sm t t- Therefore, 
one can parameterise a possible beyond the SM contribution by a quantity 
A', defined as 

Br l = Br l \ SM (l + A l ). (7) 

In the lowest order of SM the leptonic decay width of the r is due to the 
tree level W exchange, see Fig. El (left). Including the W-propagator effect 
and QED radiative corrections, the following results for the branching ratios 
in the SM are obtained (see also Sect. 4): 

Br e \ SM = (17.80 ± 0.07)%, Br^\ SM = (17.32 ± 0.07)%. (8) 

Together with the experimental data this leads to the following estimations 
for the possible beyond SM contributions to the considered branching ratios, 

A e = (0.20 ±0.51)%, A" = (0.26 ± 0.52)%. (9) 

Using them we derive the 95% C.L. bounds on A', for the electron and muon 
decay mode, respectively: 

(-0.80 < A e < 1.21)%, (-0.76 <A^< 1.27)%. (10) 

One can see that the negative contributions are constrained more strongly 
that the positive ones. 
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4 Leptonic r decays in 2HDM 



In the SM the leptonic r decay, r — > WiV T , proceeds at tree-level via the 
W ± exchange. The formula below describes this contribution in the Fermi 
approximation, with leading order corrections to the W propagator, and 
dominant QED one-loop contributions, 



tree 



G 2 F m 5 



mf_ 

2 



lf(— L) ( 1 
25 



192vr 3 m 
a\m 

1 + 



Ilr 
2tt V 4 



3m 2 
hm 2 w 



71 



rat 
2^- 



ni 



w 



x 



(11) 



f(x) = l-8x + 8x d -x 4 -12x 2 lnx. (12) 

We will denote the SM contribution in short as T l , skipping here and 
below the superscript I if not necessary. 





Figure 3: Tree-level contributions to the r leptonic decays. The W exchange 
in the SM (on the left) and the exchange in 2HDM (on the right). 



In 2HDM there is, in addition, a tree contribution due to the exchange of 
the charged Higgs boson, Fig. El (right). This new contribution is given by 



where 



tree 



m^mf tan 4 (3 



4M 4 ± 



rriim T tan 2 3 mi 
2 —r-, K 



M 2 H± 



m. 



mf 



mi 



(13) 



k(x) 



9{x) 



g(x) = 1 + 9x - 9x 2 



+ 6x(l+x)ln(a;). (14) 



Note that the second term is coming from the interference with the SM 
amplitude and it is much more important than the first one, that is sup- 
pressed by a factor m 2 tan 2 3/8Mj I± . Note that such suppression can be 
compensated by a very large tan 3 only. 

In 2HDM there are also one-loop contributions involving neutral as well 
as charged Higgs and Goldstone bosons. All these contributions are included 
in the Gp scheme as follows: 



S 



r'(i + 5Z Lt + 5Z Ll + 5Z LUt + 5Z Lui ) + r 





(15) 



where the first term corresponds to the SM prediction, Z L f = 1 + SZ L f are 
the renormalisation constants for the left component of the fermion / and 
Ffoops corresponds to the one-loop corrections to the W ± exchange tree-level 
amplitude. The H ± exchange tree-level contribution and its one-loop and 
counterterm corrections are described by r^ ee , Tf Q ^ ps and Tq T , respectively. 

The tree-level H ± contribution is numerically small and the radiative 
corrections to this amplitude will be neglected here. Taking this into account 
we will just consider the tree-level contribution eq. ()13)) . implying that 



5 One-loop 2HDM(II) corrections 

We evaluate, in the 't Hooft-Feynman's gauge, the one-loop contributions 
coming from the 2HDM(II) to the quantities A', using definitions and con- 
ventions for one-loop integrals of [21]. We will take into account the fact that 
the E ± and W ± masses are very large compared with the leptonic masses 
and external momenta, and we will neglect masses of muon and electron in 
the loop calculation. This means that the obtained one-loop corrections are 
universal, i.e. they do not depend whether decay into e or /i is considered, 
so A oneioo P = ^Ineioop = ^oneioop- Moreover, we will focus on large tan/? 
enhanced contributions. 

5.1 Renormalisation constants 

In order to evaluate the 2HDM contributions to the fermion fields renor- 
malisation constants, one has to compute the self-energies coming from the 
diagrams shown in Fig. 0J 



Figure 4: Two-point diagrams contributing to the fermion fields renormali- 
sation. Here x° = K H, A, G° and x + = H + , G + . 



r 



loops — 1 CT — 



0. 



(16) 
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Charged lepton self-energies. There are two kinds of contributions, one 
involving the exchange of a neutral boson and a second one involving a 
charged boson. The latter ones are numerically negligible since they are 
proportional to mf/M^y and mf/M] I± (for x + — G + and H + , respectively). 
Therefore we will consider the corrections coming from neutral Higgs and 
Goldstone bosons only. Since these corrections are proportional to m 2 we 
will take into account just the contributions to the self-energy of r. We 
obtain 

5Z Le = 5Z Lll = 

5Z Lr = A^ + Af + A^ + Af 

Ah G F ml sin 2 a „. 9 , r9 9x 

-H- Gprnl cos 2 a „, 9 , 9 . 

A^ = -^tan 2 /3SK;M 2 ,m 2 ) 

Af = -^S(m 2 ;Ml,m 2 )^0, (17) 

where we use the following abbreviation 

BQ = [B + B 1 + Am 2 T B' + 2m 2 T B[](). 
The G° contribution will be neglected since it is not tan 2 f3 enhanced. 

Neutrino self-energies. In this case only the H + and G + contributions 
are involved and, since again these corrections are proportional to the mass 
of the lepton in the loop, we will just consider the corrections to the tauonic 
neutrino field renormalisation. We obtain 



SZ LUe = 5Z Lv ^ = 

5Z Lvr = Af +Af 

Af = -^tan 2 (3[B + B 1 mM^,ml) 

Af = ^[Bo + B^M^mD^O. 



5.2 One-loop three-point contribution 

The one-loop three-point diagrams contributing to A in the 2HDM(II) are 
presented in Fig. We use here the following notation: x° = K H, A, G°, 

X + = H+, G+ and (V, 0) = (G+, Z), (W + , h), (W + , H)/(Z, G+) 
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T 




Figure 5: Three-point diagrams contributing to the W ± ti , t vertex correction. 
Similar diagrams exist for the W lu\ vertex. x° — h, H, A, G°, x + — H + , G + 
and (V, <P) = (G+, Z), (W + , h), (W + , h)/(Z, G+). 



These W lv\ vertex corrections are proportional to the lepton mass and 
therefore we will consider only the radiative contributions to the W ± tv t 
vertex. The different contributions coming from each diagram are as follows. 



X + — X° ~ T Loops. We have computed them (Fig. left) in the limit of 
large Mjj± and Mw That means, we have obtained the complete expressions 
and keep only such terms that do not decouple in the limit M H ±, Mw 3> m T . 
The resulting expressions are: 



H+h 
loops 

H+A 
loops 

H+H 
loops 

A? +h 



A 



A 



A 



G F m 2 r 
G F m 2 T 
G F m 2 T 

2y/2<K 2 



tan/5 -cos(a - i3)C 20 (m T , m v ; M h} m T} M„±) + 

cos/3 t ' 

tan 2 (3C 20 (m 2 T , m^; M\, m 2 , M 2 H± ) + ... 



cos a 



tan/3 -sin(a - (3)C 20 (m 2 , mi^Mfj, m 2 , M#±) + ... 



a g+h 

loops 



/\G+G° 
loops 



Gpm 1 . sin a; 
2y/2n 2 cos/3 
Gftti 2 cos a 
2V2-K 2 cos/3 
G F m 2 



2V2-K 



cos P 

sin(a - /3)C 20 (m 2 , m 2 r ; M 2 , m 2 , M^) + . . . ~ 
cos(a - p)C 2 o(m 2 T , m 2 r ; M 2 H , m 2 , M w ) + . . . ~ 

(19) 



r -C 20 (m 2 T , ml; M 2 ZJ m 2 T , M 2 W ) + . . . ~ 



The three last contributions can be neglected in the large tan/3 limit. 



V — <fi — I Loops. These contributions (Fig. 03 middle) are numerically 
negligible as do not contain any tan 2 f3 factor. Therefore, in our work 



A 



V(j> 
loops 



0. 



(20) 
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Figure 6: Dominant one-loop diagrams contributing to leptonic r decays in 
the limit of large tan (3, here 4>° = h, H, A. 



t — v T — x + Loops. We have computed these contributions (Fig.EJ right) 
and checked that they decouple in the limit of very heavy charged Higgs 
boson and boson, as their leading terms in this limit are proportional to 



or Therefore 



0. (21) 



5.3 One-loop box diagrams 

The one-loop box diagrams also contribute to the r leptonic decays. All 
of these diagrams involve the exchange of a charged Higgs boson or a 
boson. They can be safely neglected due to the mass dimension of the D 
integrals that describe these diagrams, namely 

D ~ — r , D u ~ , D, w ~ — , D uvo ~ — , D wcri ~ O(M ). (22) 

Since Mh± and M^/ are very large as compared to m T , they will drive 
the mass dependence of the integrals, so M = Mh±, Myy. Therefore only the 
terms proportional to D^ vpi do not decouple and give relevant contributions. 
However, in the considered case of r decays there are no such contributions. 
Therefore we can neglect box diagrams altogether: 

At * (23) 



5.4 Final expression for one-loop contribution 

Taking all this into account, the dominant diagrams in the limit of large 
tan j3 are reduced to the ones drawn in Fig. |U1 The contributions coming 
from these diagrams are 



G F m 2 T 2 

^oneloop = g^/^ 2 ^ n $ * 
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[- cos 2 (P - a)B(m 2 T ; M 2 , m 2 ) 

-B(m 2 T ;M 2 A ,m 2 T ) 

-sin 2 (/3 -a)B(m 2 T -M 2 H ,m 2 T ) 

-2[B + B 1 ](0;M 2 I± ,m 2 r ) 

+4 cos 2 (/3 - a)C 20 (m 2 T , ml T ;M 2 h , m 2 , M 2 H± ) 

+4C 20 (m 2 ,m 2 T ;M 2 ,m 2 ,M| ± ) 

+4sin 2 (/3 -a)C 20 (m 2 T ,ml T -M 2 H ,m 2 T ,M 2 H± )} . (24) 

An easy to handle expression can be obtained from eq. ()24j) for neutral 
Higgs masses larger that the r mass, M^o > m T . Notice that no assumption 
on the Higgs spectrum is made 5 . In this limit, we get 6 



^oneloop 



G F m 2 T 
8V2tt 2 

A = 



tan 2 /? A 



+- 



In 
In 



m 2 . 
Ml 
m 2 



+ F(R H ± 

+ F(R A 
(Ml 



+ -cos 2 (/?-«) (ln( )+F(R h ] 



+ -sin 2 (/3 - a) ( In 



M| 
m 2 



+ F(R H ] 



(25) 



where R^ = M^/M^ and 



F(R) = -1 + 2 



R 2 \nR 2 



1 -R 2 

Some useful limits of the F function are: 



F(R<£1) 



F(R 



-3 , F{R > 1) 



1 + 21niT 



(26) 



(27) 



The above expression depends logarithmically on the ratios of the mass of 
each Higgs boson to the mass of tau lepton. However, this may be misleading 
as in fact there is no dependence on mass of the tau lepton in the above 
expression. Indeed, one can see this by looking at other useful form of A, 7 



3 + - (G(R A ) + cos 2 ((3 - a)G{R h ) + sin 2 (/? - a)G(R H )) 



(2f 



5 We generalised here the result in [9] where assumptions were made M#± , Ma 3> Mh 
and a — (3. 

6 In agreement with [22] result derived in the context of the MSSM. 
7 We thank M. Misiak for this suggestion. 
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where 



G(R) =\nR 2 + F(R). 



(29) 



In the following, to explore the phenomenological consequences of the 
large tan f3 enhanced 2HDM(II) radiative corrections to the leptonic r decays, 
we will use both the exact and approximated expressions (|2^|) and (|25H28|) . 
respectively. 

5.5 One-loop corrections for some interesting scenarios 

In some phenomenologically interesting scenarios, the expressions (J25I28J) can 
be simplified. In the case of light h and sin 2 (/3 — a) = 0, A does not depend 
on Mh and two limits are worth to be studied, 

M A = M H ± — > A = ln^^+1 and M A < M H ± -> A = ln^^+ln^-+2. 

(30) 

Notice, that when h does not couple to gauge bosons and therefore the Hig- 
gsstrahlung process at LEP is not sensitive to such Higgs boson, the leptonic 
tau decays have maximal sensitivity to h as A depends logarithmically on 
its mass, without any suppression factor. 

If A is light and sin 2 (/3 — a) = 1, one obtains the same expression for A 
that in the previous case with obvious replacing h — > A and A — > H . There- 
fore any analysis with h light and sin 2 (/? — a) = can be easily translated 
to the case of light A and sin 2 (/5 — a) = 1. 

The useful expression which holds for arbitrary sin(/3 — a) and degenerate 
H, A, H ± (with a common mass M) is: 

A = cos 2 (/3-a) [ln^ + 1]- (31) 

We see that in a SM-like scenario, with light h, sin 2 (/? — a) = 1 and very 
heavy degenerate additional Higgs bosons, A goes to zero, what signals a 
clear decoupling. 

6 Numerical analysis 

In this section we analyse the dependence of the 2HDM(II) one- loop cor- 
rections obtained in the previous sections for the leptonic r decays on the 
different Higgs bosons masses and mixing angles. First we stress that typi- 
cally the one-loop contribution dominates the 2HDM(II) effects. They are, 
for fixed value of large tan f3 and in the interesting region of parameter space, 
five orders of magnitude larger than the corresponding tree-level H ± contri- 
bution to A e , and one or two orders of magnitude larger for the A M . There- 
fore, although we will include all contributions in the numerical analysis, the 
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main features of the 2HDM effects are described by the one-loop correction 
in eq. (J25|) . In the following only results for muon decay channel will be pre- 
sented. We stress however that the obtained one-loop corrections are the 
same for the electron and muon channels. 

In eq.()25|) one can distinguish two contributions, one coming from the 
charged Higgs alone and the other one involving also the neutral Higgs 
bosons. The former is always negative and it becomes more negative for 
a larger charged Higgs mass. The latter is typically positive and it grows 
with the neutral Higgs masses. In this way the total 2HDM(II) one-loop 
effects, being a sum of two contributions of the same order and with different 
signs, will be large when one of these contributions dominates. Since the 
modulus of both corrections grow with the Higgs masses one expects large 
one-loop effects in two cases: (i) heavy if ± and light <p° (large negative cor- 
rections) and (ii) light and heavy <p° (large positive corrections). Taking 
into account the lower bound for Mh± coming from b — *> S7, above 490 
GeV, one expects to get large radiative effects in case (i) only. Note that in 
this (i) case the A (loop contribution) is negative, as well as the tree-level 
H ± exchange, eq. tflHj) . 

We will focus on two scenarios of a special phenomenological interest: 
with a light scalar h and with a light pseudoscalar A. As all contributions 
considered here are proportional to tan 2 (3, they will be plotted for tan [3=1, 
to be rescaled by tan 2 (3. 

6.1 Light scalar Higgs boson, h 

First we will consider a scenario with a light scalar boson, h, with mass 
Mh below 114 GeV, and degenerate heavy Higgs bosons, with masses Ma = 
Mh = Mh+ = M, above 300 GeV. For such a light Higgs boson h, its 
couplings to gauge bosons are constrained by LEP data as shown in Fig. Q 
(left), lying between and sin 2 (/3 — a)\ max . Note that for arbitrary sin(/9 — a) 
and degenerate H, A, H ± the equation (}3*T]) holds. 

In the light h scenario instead of the degenerated heavy additional Higgs 
bosons, one can also consider a spectrum with SM-like Higgs boson H ( i.e. 
with couplings to the gauge bosons as for the SM Higgs, namely Xv = 1 or 
sm(/3 — a) = 0). It is reasonable to assume that such Higgs boson has a mass 
in the region expected for the SM Higgs boson, say Mh = 115 GeV, although 
as follows from eq. f)25j) nothing depends on this mass. On contrary, one gets 
here a clear dependence on mass of h: 

A = ln^ + 1. (32) 

The different contributions to A oc A tan 2 /? are plotted in Fig. (left) 
for Mh = 5 and 70 GeV, for degenerate heavy Higgs bosons. The total 
( ie.sum of the tree and one-loop) contributions are plotted using solid 
lines, while the one-loop contributions using dashed lines, respectively. As 
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can be seen, the H ± tree-level effect is important for low M but the one- 
loop contribution becomes dominant for M > 500 GeV. In particular, the 
logarithmic dependence on M coming from the one-loop corrections is clearly 
seen. Notice that curves are plotted for sin 2 (/3 — a) = and sin 2 (/? — a) \ max , 
the maximum value allowed by LEP data for a given Mh value. For h mass 



equal to 5 GeV the results for different sin 2 (/3 
0.02, can not be distinguished. 
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Figure 7: The total (solid line) and one-loop (dashed line) contributions to A 
for tan /3 = 1. Left: Results for Mh = 5 and 70 GeV are plotted as a function 
of M = M A = M H = M H ± . Results for mass of 70 GeV for sin 2 (/3 - a) = 
and 0.3 are presented by the bottom and upper lines, respectively. Right: 
Results for M A = 100 GeV, M H ± = 4 TeV and sin 2 (/5 - a) = (M H is 
arbitrary), as a function of Mh are shown. 



The dependence of A on the light Higgs mass can be seen in Fig. [7| (left) 
by comparing the results obtained for Mh = 5 and 70 GeV. This dependence 
is explicitly presented in Fig. [7| (right) where the contributions are plotted 
as a function of M h , for M A = 100 GeV, M H ± = 4 TeV, sin(/5 - a) = 0. 
The 2HDM(II) one-loop corrections decrease logarithmically with increasing 
Mh as described by eq. (|52~j) . So the lighter h the larger the one-loop correc- 
tions. One can see that A decreases linearly with increasing sin 2 (/3 — a), in 
agreement with eq.(pTT|). 

In the case with the SM-like H we have sin 2 (/3 — a) = 0, then A becomes 
insensitive to the value of M H , see eq.(|25j) and discussion above. Therefore 
all above results obtained for the sin 2 (/3 — a) = case hold also here for the 
SM-like H. 
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6.2 Light pseudoscalar, A 



In the case where the pseudoscalar, A, is light and the two neutral scalars 
are degenerate Mh = Mh, the A does not depend on sin 2 (/? — a). For 
Mh = Mh = Mh± = M we get a simple formula 

A = ln^ + 1. (33) 

It is similar to the formula obtained for the discussed above case of light h 
for sin 2 (/? — a) = 0, with obvious change Mh to Ma- Therefore we will not 
present the results corresponding to such light A case. 

There is an interesting light A scenario where in addition to A also h is 
not very heavy. We call this case a light A & h scenario. Here we choose 
the h mass to be equal to 100 GeV to avoid a direct conflict with LEP 
data presented in Fig. Q (right). In Fig. [H] the total contribution to A is 
plotted as solid lines, while one-loop corrections as dashed lines, respectively. 
Also in this light A & h scenario we see that the one-loop effects dominate 
for large M scale. The largest deviation from the SM prediction occurs for 
sin 2 (/? - a) = 0. 

Contribution to Br_mu in 2HDM(II) 
-5e-07 — 1 




-3e-06 1 1 — 1 — — 

1000 

M, GeV 



Figure 8: The total (solid line) and one-loop (dashed line) contributions to 
A for tan (3 = 1. Results for M A = 40 GeV and M h = 100 GeV are shown 
for sin 2 (/3 — a) = and 1, bottom and upper lines, respectively. 



6.3 Comparison of the exact and approximated results 

Results based on eq.(J2H) and the approximation eq. (|2H|) have been plotted 
together in all the figures, being clearly indistinguishable. Therefore, the sim- 
ple approximated formulae eqs. (|25l28j) can be used to describe the 2HDM(II) 
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one-loop corrections to the r leptonic decays in the whole considered range 
of parameters. 

7 Constraining 2HDM(II) by the r decay data 

In this section we use the leptonic r decay data to constrain 2HDM(II) 
parameters. The complementarity between LEP processes used for direct 
searches of light Higgs bosons and indirect measurements based on the lep- 
tonic r decays will be exploited to explore the "pessimistic", for the direct 
searches, scenarios. In particular the case sin 2 (/3 — a) = will be studied, 
since in this scenario the Higgsstrahlung and VV fusion processes for h are 
suppressed. The 95 % CL bounds for A derived by us in Sect. 3 allow to set 
upper bounds on tan j3 (Yukawa couplings) for both the light h or A scenar- 
ios. We provide also exclusion for (M^, Ma) plane for various values of tan/? 
and sin(/3 — a). 

In addition, we obtain from r decays new lower and upper bounds on the 
charged Higgs boson mass as a function of tan (3. 

7.1 Constraints on the Yukawa couplings of the light- 
est neutral Higgs bosons 

The upper limits on tan/3 (Yukawa coupling Xd) for light h and light A 
scenarios are shown in Fig. El and EH respectively. 

In the "pessimistic" light h scenario with sin 2 (/3 — a) =0, the leptonic 
r decay data can be exploited to set upper limits on the Yukawa couplings 
as a function of M^. They can be compared with limits coming from other 
experiments. 

In Fig. El the upper limits on the tan /3 (Yukawa couplings) for light 
h, assuming sin(/3 — a) = 0, derived from the leptonic r decay data are 
presented. One can see that these data provide upper limits on tan/3 in 
region unaccessible by other experiments, namely for mass above 45 GeV. 

As an opposite case to the light h scenario, one can consider the case 
with a light pseudoscalar A. If sin(/3 — a) = 0, M4 can be low if h is heavy 
enough to suppress the associated (h, A) production. The Yukawa couplings 
of A can be then constrained just by the Yukawa process with ffA final state 
and the Upsilon decay, T — > Aj, for a very light Higgs boson A. Also in this 
case the leptonic r decays can be used to set upper limits on the Yukawa 
coupling (tan/3) as a function of Ma, see Fig. [TUl The right panel shows the 
region around mass of A equal 10 GeV. 

Since this scenario can be relevant in explaining (g — 2) M data, we plot in 
Fig. ^2 the upper limits for tan/3 from the leptonic r decay and the allowed 
region from the newest g — 2 for muon data, and for a comparison all other 
existing upper limits for A. Degenerate masses of h, H, H +s were assumed 

8 sin(/3 — a) is then arbitrary 



18 



Exclusion 95%C.L. for h in 2HDM(II) 
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Figure 9: 95 %CL upper limits from r decay for tan/3 and sin(/3 — a)=0, 
as a function of Mh compared to the existing upper limits from the Yukawa 
processes at LEP (Opal, Delphi) and the Upsilon decay. The two almost 
horizontal lines (in green) corresponds to Ma = 100 GeV, for Mh± = 500 
GeV and 4 TeV, upper and lower lines, respectively. The results for the 
degenerate A and H + with mass 4 TeV are plotted by using thicker line. 




to be equal to 1 and 4 TeV, the corresponding upper limits are presented by 
the upper and lower lines, respectively. 

7.2 Constraints on a light A & h scenario 

A scenario with both h and A light is also of phenomeno logical interest. Since 
A can be large for low Mh and Ma, the leptonic r decay data can be used to 
constrain this scenario, in the (Mh, Ma) parameter space. The comparison 
of these constraints with the ones coming from direct searches will reveal the 
importance of indirect ones from the leptonic r decays. 

In Fig. [^1 the constrained regions in the (Mh, Ma) plane, laying between 
axes and the corresponding curves, are shown for sin(/3 — a) = and tan f3 
equal 60 and 90. The excluded regions are symmetric in Mh and Ma] they 
rule out the possibility of both h and A being very light. These constraints 
should be compared to the constraints shown in Fig. ^ (right). For large 
values of M^± and tan/3 the 2HDM(II) one-loop effects can be very large 
and some of the regions of the parameter space allowed by direct searches 
can be excluded indirectly by using the leptonic r decays. 
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Figure 10: 95 %CL upper limits from r decay for tan (3 as a function of 
compared to the existing upper limits from Yukawa processes at LEP ( OPAL, 
DELPHI) and Upsilon decay. The two almost horizontal lines (in green) 
corresponds to Mh = 100 GeV and M H ± = 500 GeV and 4 TeV, upper and 
lower lines, respectively. The results for degenerate h and H ± with mass 4 
TeV are plotted in thicker line. Left: Mass range for A from 5 to 200 GeV, 
Right: Mass range for A from 8 to 12 GeV. 



7.3 Constraints on the charged Higgs-boson mass 

From the leptonic tau decays one can derive limits on the mass of charged 
Higgs boson as a function of tan (3. A standard derivation within the 2HDM 
(II) is based on the tree-level H + contribution for the leptonic tau decay into 
a muon. Such derivation can be found in almost all papers, both theoretical 
and experimental ones, devoted to this subject (see eg. [1]). 

First, we apply such a standard method to derive from the tree-level 
contribution (for muon) only the lower mass limit for H + . By applying 
the obtained lowest value for the 95% CL deviation from the SM prediction 
(eq. [TUJ) we updated the existing lower mass limit. We got the following limit 

M H ± > 1.71 tan/? GeV (34) 

with coefficient 1.71 to be compared to the corresponding coefficients from 
[23] and [24], equal to 1.86 and 1.4, respectively. Note, that this is nothing 
else, up to the lepton mass ratio, what issued as the constraints on the Michel 
parameter 7] in the 2HDM (II), see eg. [1,25] 

Next, knowing that the one- loop corrections are typically more impor- 
tant than the tree contribution, we use them in the derivation of the mass 
limit for H + . We observe that since A oneloop grows with M H ±, this one-loop 
correction allows to put upper bounds on M^± in scenarios with light neutral 
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Exclusion 95%C.L for A in 2HDM(II) 
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Figure 11: Upper limits for tan/? from the leptonic r decay (thick grey lines) 
and the allowed region from the newest g — 2 for muon data (thick black 
lines), in comparison all other existing upper limits as a function of Ma- 
Degenerate masses of h, H, H + were assumed to be equal to 1 and 4 TeV; the 
corresponding results from tau decay are given by the upper and lower thick 
grey lines, respectively. 

Higgs bosons. In particular, for sin(/3 — a) = 0, A goes as \n(M H ±/Mh) + 
ln(M#± / 'Ma) (eq. EUJ) and therefore the lighter h and A the stronger upper 
bounds for M H ±. 

In Fig. Uni (left) the individual lower and upper bounds on M H ±, as ob- 
tained from the tree-level if ± exchange diagram only and from the one-loop 
contribution only, are plotted as a function of tan j3. The constraints based 
on the one-loop contributions are plotted for various masses of h, equal to 
5, 20 and 100 GeV, assuming sin(/3 — a) = and a degeneracy in masses of 
A and H^. Upper limits coming from the one-loop corrections are plotted 
both for the muon and electron decay channels, the limits from electron one 
is slightly weaker (dashed (green) lines). The relevant lower limits from the 
tree-level H ± contribution is obtained only from the muon channel. The 
lower bound coming from b — > S7 analysis is also shown for a comparison. 

In Fig. H31 (right) we present results as described above for one particular 
mass of h equal 20 GeV, together with a full bound based on a sum of tree 
and one-loop contribution (thick line). We see that a full bound gives both 
the lower and upper limits as a function of tan (3. In this figure we present 
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95 % CL exclusion Mh & MA region 




20 40 60 80 100 120 
Mh, GeV 



Figure 12: The excluded regions in the (M h , Ma) plane for sin 2 (/? — a) =0. 
The excluded regions lay between axes and curves corresponding to Mh± = 
500 GeV and 1 TeV, dashed and solid lines, and to tan/3 = 60, 90, down and 
upper lines, respectively. 



also results obtained for other mass of A, equal 100 GeV. The thin black 
dotted line corresponds to the limits obtained then, with all other 2HDM 
parameters as used above to obtain a thick line. 

It is clear that the constraints for mass of H + change drastically if the 
one-loop contributions are included in the analysis. In particular, the lower 
bound is higher than the tree-level limit, eq.([34|). Only for the SM-like h 
scenario, with sin(/3 — a) = 1 and all other Higgs boson mass heavy and 
degenerate, the tree-level contribution gives a reliable estimation. 

Based on results presented in Fig. EH the restrictions can be set on M#± 
for large values of tan/3 (tan/3 > 60). In particular, in a scenario with light 
h and not so heavy A, M#± should be lower than 3 TeV for tan/3 = 65. 
Although large values of tan (3 are required, this upper bound to the charged 
Higgs mass is important due to the difficulty on setting upper bounds on 
masses of undiscovered particles. 

8 Summary and conclusions 

In this work we have computed the 2HDM(II) one-loop corrections to the 
leptonic r decays. As a main result we have obtained that these one-loop 
effects are larger than the corresponding tree-level contribution in the 
relevant regions of the parameter space. Our analysis has been focused on 
the tan/3 enhanced contributions and an easy-to-handle formula has been 
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Figure 13: Limits on charged Higgs boson mass as a function of tan (3 ob- 
tained from leptonic tau decays. The upper limits from the one-loop contri- 
bution (dashed and solid lines correspond to the electron and muon channels) 
and lower limits (straight lines) from the tree H + exchange from muon chan- 
nel are shown. Lower limit from b — * 57 is also shown. Left: The upper 
limits obtained from the one-loop corrections only for Mh = 5,20,100 GeV 
and sin 2 {(3 — a) = 0, assuming Ma = M^ are presented. Right: The same as 
in Left for one mass = 20 GeV, in addition in form of the full constraint 
from the total (loop plus tree) contribution, for degenerate masses Ma = M H ± 
(thick solid line), and for Ma = 100 GeV (thin dotted line) is presented. 



obtained describing these one-loop effects in the approximation of the Higgs 
boson masses larger than the r mass. This formula allows us to study all the 
2HDM parameter space in a transparent way. 

After the numerical analysis of the corrections, the constraints on the 
2HDM(II) parameters from the leptonic r decay data have been obtained 
in different scenarios. In particular the "pessimistic" scenarios for the direct 
searches of light h Higgs bosons at LEP have been intensively analysed. From 
this analysis we have obtained upper limits on the Yukawa couplings for both 
light h and light A scenarios, constraining also the light Ahh scenario. We 
have updated the existing in literature lower limits on Mh±, different from 
the ones coming from tree-level exchange only, and we have also obtained 
interesting upper limits on M#± as a function of tan (3. 

Therefore, one can conclude that leptonic r decay data constrain 2HDM(II) 
scenarios with large tan /?, heavy and light neutral Higgs bosons. 

Obviously, the large 2HDM(II) one-loop corrections found in this paper 
can have consequences for other type of processes, which will be analysed 
elsewhere. 
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